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ABSTRACT 

We present forecasts for constraints on deviations from Gaussian distribution of primordial 
density perturbations from future high-sensitivity X-ray surveys of galaxy clusters. Our anal- 
ysis is based on computing the Fisher-Matrix for number counts and large-scale power spec- 
trum of clusters. The surveys that we consider have high-sensitivity and wide-area to detect 
about 2.5 x 10 5 extended sources, and to provide reliable measurements of robust mass prox- 
ies for about 2 x 10 4 clusters. Based on the so-called self-calibration approach, and including 
Planck priors in our analysis, we constrain at once nine cosmological parameters and four nui- 
sance parameters, which define the relation between cluster mass and X-ray flux. Because of 
the scale dependence of large-scale bias induced by local-shape non-Gaussianity, we find that 
the power spectrum provides strong constraints on the non-Gaussianity /nl parameter, which 
complement the stringent constraints on the power spectrum normalization, cr 8 , from the num- 
ber counts. To quantify the joint constraints on the two parameters, era and /nl, that specify 
the timing of structure formation for a fixed background expansion, we define the figure-of- 
merit FoMsft = (det \Cov(cr%, ,/nl)]) - • We find that our surveys constrain deviations from 
Gaussianity with a precision of A/nl — 10 at lcr confidence level, with FoMsft - 39. We 
point out that constraints on /nl are weakly sensitive to the uncertainties in the knowledge of 
the nuisance parameters. As an application of non-Gaussian constraints from available data, 
we analyse the impact of positive skewness on the occurrence of XMMU-J2235, a massive 
distant cluster recently discovered at z — 1 .4. We confirm that in a WMAP-7 Gaussian ACDM 
cosmology, within the survey volume, 5 x 10~ 3 objects like this are expected to be found. 
To increase the probability of finding such a cluster by a factor of at least 10, one needs to 
evade either the available constraints on /nl or on the power spectrum normalization cr 8 . 



1 INTRODUCTION 



The standard inflationary scenario, based on the single scalar 
field slow-roll paradigm, predicts primordial density perturbations 
to be virtually indistinguishable from a Gaussian distribution. 
However, a number of variants of inflation have been proposed 
which are able to generate a certain amount of non-Gaussianity 
(e.g.. iBartolo et alj|2004 ICherj|20ld) . Therefore, testing to what 
precision we can measure possible deviations from Gaussianity 
with available and future observations has important implications 
on our understanding of the mechanism that seeded density fluctu- 
ations in the early Universe. Analyses of the Cosmic Microwave 



Background (CMB) provide at present the tighest constraints on 
the amount of allowed non-Gaussianity. A number of analyses 
based on the WMAP data conv erge to indicate consistency with 
the Gaussian ass umption (e. g., Komatsu et al. 2010 and references 
therein; cf. also [Yadav & WandelJ l2008h . While data from the 
Planc k satellite are expected to further tig hten such constraints 
(e.g., lYadav et al.l 120071 : iLiguori et al.l 120101) . it is worth under- 
standing whether non-Gau ssianity can be probed by larg e-scale 
structure observations (e.g., Slosar et a lj2008l : IVerd e2010). 



Non-Gaussian perturbations are expected to leave their im- 
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print also on the pattern of structure growth at least in two different 
ways. First, we expect that a positively skewed distribution pro- 
vides an enhanced probability of finding large overdensities. This 
translates into an enhanced probability of forming large collapsed 
structures at high redshift, thereby changing the timing of struc- 
ture formation and the shape and evolution of the mass function of 
dark-matter halos. 

After the first pioneering studies of the effect of non- 
Gaussiantiy on the mass function (e.g. | Matarrese et alj 1 19861 : 
IColafrancesco et alj fl98^ : H organi & Bo nomettol ll990). a num- 
ber of anal yses have been carried out since the beginning of 
2000s (e.g.. IMatarrese et alJl200Ct iMathis et alj|2004 iKang et al l 



2007 ; ISefusatti et ai]l2007l : iGrossi et al]|2007l ; iMaggiore & Riottol 
2009). Furthermore, the realization of large-scale cosmologi- 



cal simulations with non-Gaussian initial conditions have re- 
cently provided a validation of the non-Gaussian correction to 



20071: iDalal et al.ll2008l:lDesi acmes et alj|2009: 


Grossi et alj|2009l: 


Giannantonio & Porciani 20091: Pilleoich et al. 


201Ch. More re- 



cently non-Gaussianity effects on the large-scale distrib ution of col 



lapse d halos were studied: it has b e en demost r ated ( e.g., Dalai et al 



2008 1: IMatarrese & Verdd 120081 : IValageasI l20ld : lLam & Shetti 
2009) that non-Gaussianity affects the large-scale clustering of ha- 
los in such a way that the linear biasing parameter acquires a scale 
dependence. This modifies, in a detectable way, the power spec- 
trum of the distribution of any tracer of cosmic structures at small 
wavenumbers and offers a unique way of testing the nature of pri- 
mordial fluctuations. 

The evolution of the mass function of galaxy clusters identi- 
fied in X-ray surveys has been e xtensively used in the past to con 



strain cosmological models (e.g., Borgani et alj|200ll : iRosati et al.1 



120021 : [Schuecker ~e7aT1l2003l : Voit 2005). These studies have re 
cently attracted renewed interest, thanks to detailed follow-up ob- 
servations of clusters selected from ROSAT observations, they 
have been carried out either to constrain the Dar k Energy equa- 
tion of state within the G aussian paradigm (e.g., Ivikhlinin et al.1 
2009; iMantz et alj 120091), or to tes t poss i ble deviations from 
standard gravity (e.g.. ISchmidt et al.1 120091 : iRapetti et all l2009h . 
Ijimenez & Verdd J2009I) have recently analyzed the effect of non- 
Gaussianity on the population of massive high- redshift clusters, 
like the one recently discovered by Jee et al. (2009) at z — 1.4. It is 
important to remember that the cosmological constraints obtained 
from clusters so far have been derived from small ROSAT-based 
samples, containing ~ 100 clusters at z < 1. It is therefore easy to 
imagine the vast margin of improvements offered by next genera- 
tion X-ray surveys which will detect ~ 10 5 clusters to z ~ 2. 

These surveys should cover a large enough volume at high 
redshift to test non-Gaussianity in the regimes where its effects 
are clearer, namely the high-mass tail of the mass fun ction and the 
large- scale power spectrum o f the cluster distribution. iFedeli et"al] 
( 2009) and Roncarelli et alj d2010h presented predictions for the 
number counts and clustering of galaxy clusters expected from the 
eROSITA X-ray survey and from the Sunyaev-Zeldovich SPT sur- 
vey. While these analyses confirmed the potential of these surveys 
to provide interesting constraints on non-Gaussian models, they did 
not include detailed forecasts on the constraints on non-Gaussian 
models and a detailed assessment of the effect of uncertainties 
in the scalin g relations between cluster masses and observables. 
lOguri ( 2009) foll owed the self-calibration approach by Lima & HvJ 
(2005) (see also lMaiumdar & Mobj|2004l : lBattve & WellefeOoH) 
to forecast the capability of future optical cluster surveys to con- 
strain non-Gaussian models. This study showed that combining 



number counts and clustering of galaxy clusters can potentially pro- 
vide quite strong constraints on deviations from Gaussianity. 

The aim of this paper is to derive forecasts, based on the 
Fisher-Matrix approach, on the capability of future X-ray clus- 
ter surveys to constrain deviations from Gaussian perturbations. 
Besides focusing on the characteristics of X-ray, rather th an op- 
tically selected samples, our analysis differs from that by Oguri 
(2009) for the meth od to include information fr om large-sca l e clus- 
tering. lOguril j2009h adopted the approach by Lima & Hu d2005t) 
where clustering is included by accounting for fluctuations of 
cluster counts within cells having a fixed angular size. This im- 
plies that, at each redshift, clustering information is restricted to 
one physical scale. In our analysis , we follo w the approach origi- 
nally presented by Tegmarkl J 19971) (see also lFeldman et alii 19941 : 
iMaiumdar & Mohj 120031 7 in which the clustering Fisher-Matrix 
is computed for the allowed range of wavenumbers, by weighting 
them according to the effective volume covered by the surveys. 

Another distinctive aspect of our analysis is that it is based 
on X-ray surveys of next generation, whose sensitivity and angular 
resolution are high enough to warrant an accurate measurement of 
robust mass proxies, related to the cluster gas mass and X-ray tem- 
perature for a large number of clusters. As we will discuss in the 
following, surveys with these characteristics can be provided by 
an already proposed X-ray telescope, which joins a large collect- 
ing area to a large field-of-vi ew and a high angular resolution over 
the en tire field of view (e.g.. lGiacconi et al1l2009l : lvikhlinin et al.1 
2009). The great advantage of having a similar survey is that there 
is no need to assume any external follow-up observation for a sub- 
set of identified clusters. Moreover, the possibility to define a flux- 
limit down to which measuring accurate mass proxies for all clus- 
ters allows one to set robust priors on the scaling relations between 
cluster mass and observables, which is one of the main source of 
uncertainty in the co smological application of galaxy clusters (e.g., 
lAlbrecht etafl 20091) . 

In principle, the method used in our analysis can be applied 
to any cluster surveys, including optical and SZ ones. Although so 
far X-ray surveys have been mostly used for cosmological applica- 
tions of clusters, upcoming large optical and SZ surveys promise 
to provide an important contribution to this field. Our method only 
requires a well defined selection function and calibrated mass prox- 
ies. Since cluster surveys at different wavelengths have different ef- 
ficiencies probing different mass ranges at different redshifts, they 
will ultimately provide complementary approaches to the deriva- 
tion of cosmological constraints. 

The plan of this paper is as follows. In Section 2 we sum- 
marize the formalism to compute non-Gaussian corrections to the 
mass function and the linear bias parameter of collapsed halos. In 
Section 3 we describe our approach to compute the Fisher Matrix 
for both the number counts and the power spectrum of galaxy clus- 
ters. In Section 4 we first describe how we compute the selection 
function and the redshift distribution expected for the X-ray sur- 
veys, then we present the results in terms of constraints on the 
parameter space defined by the non-Gaussianity parameter, / NL , 
and the power spectrum normalizations, erg. Section 5 is devoted 
to the discussion of these results. In this section we will also dis- 
cuss the competing effects of non-Gaussianity and normalization of 
the power spectrum on the expected number of clusters at z > 1 .4, 
which have a mass of, at least, 5 x 10 14 M Q , as the one recently 
studied by I Jee et"aT, I d2009h . We summarize our main conclusions 
in Section 6. 
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2 NON-GAUSSIAN INITIAL CONDITIONS 



2.1 Mass function 



Generalizations of the most standard model of inflation dGuthl 
Il98ll) give rise to seed primordial density fluctuatio ns that fol- 
low a non-Gaussi an probability distribution (e.g., iBartolo et al.l 
l2004ICherJl2010l for reviews). A particularly convenient way 
to parametrize the deviation of this distribution from the Gaus- 
sian one consists of writing the Bardeen's gauge invariant potential 
<b as the sum of a linear Gaussian term and a quadratic correc- 
tion dSalopek & Bondll99d : iGangui et al]|l994 IVerde et al.ll200Cl : 
iKomatsu & Spergell200ll) 



4> = $G + /NL * (®G - <®G>) ■ 



(1) 



In Eq. (0, the symbol * denotes convolution between functions 
and reduces to simple multiplication only in the particular case of 
constant / NL , while in general it is a function of the scale. Note 
that on scales smaller than the Hubble radius the function O equals 
minus the Newtonian peculiar gravitational potential. The funda- 
mental parameter / NL denotes the amplitude of the deviation from 
Gaussianity, and is related to the skewness of the distribution (see 
below). 



tion (as opposed to the CMB convention, see 


Afshordi & Tollevl 


2008; Pillepich et al.l20ld:ICarbone et al.l2008: 


Grossi et alj2009l) 



for defining the parameter /ml- This means, among other things, 
that the constraints given on / NL by the CMB have to be 
raised, according to the linear growth of structures, as / NL = 
g(+oo)/™ B /g(0) - 1.3/™ B , where g(z) is the linear growth sup- 
pression factor with respect to the Einstein-de Sitter cosmology. 

If the distribution of primordial density (and potential) per- 
turbations is not Gaussian, it cannot be fully described by a power 
spectrum expressed as P<t>(k) = Bk"~ 4 (where k = ||k||), but we need 
higher-order moments such as the bispectrum 6 ( i,(k 1 ,k2,k 3 ). In 
particular, different models of inflation give rise to different shapes 
of the bispectrum. In the following we adopt one particular shape, 
called local shape, means that the bispectrum is maximized for con- 
figurations in which one of the three momenta is much smaller than 
the other two ("squeezed" configurations). Inflationary models ex- 
ist that produce different sh apes for the primordial bispectrum, e.g., 



the equilateral shape (see ICrociani et all 



20091; iFedeli et al 



for applications), or t he enfolded shape Holman & Tollevl 20081 ; 



iMeerburg et alj|2009l ; IVerde & Matarresel 



2009), 



2009 



however, the lo- 

cal shape is the one giving the large st effects especially on bias 
jFedeli et alj2009l ; lTaruva et al.l2 008). hence we limit our analysis 
to this case only. 

In this case the parameter / NL is a dimensionless constant 
and the bispectrum can be written as (e.g., ICreminelli eTal]|2007l ; 
iLo Verde et alj2008h 



fi*(ki , k 2 , k 3 ) = 2/ NL 5 2 [fcT 4 C 4 + + k'r 4 kT 4 ] , 



(2) 



where n is the primordial spectral index and B is the amplitude of 
the spectrum P<i>(k), related to the amplitude A of the power spec- 
trum of density fluctuations, P{k) = Ak", by the relation 



B=-AH 4 n Cll, 



(3) 



Two of the ingredients that critically affect the observed prop- 
erties of the cluster population are influenced by non-Gaussian ini- 
tial conditions: the halo mass function and the linear bias of dark- 
matter halos. 



For the mass function of cosmic structures in non-Gaussian cos- 
mologies different prescriptions exist. The one adopted here is 
taken from lLo Verde et alj (12008). where the authors approximated 
the probability distribution for the smoothed dark-matter density 
field using the Edgeworth expansion and then performed the inte- 
gral of the probability distribution for threshold crossing exactly 
on the first few terms o f the expansion itself. The result for the 
|Press&Schechteri ( ll974l) (PS) mass function reads 



n Na (M,z) 



lit 
V ttM 



exp 



'2crl 



S 3 cr M I S 4 (z) 



+ T" 



6 

1 dS 3 



diner M I d c (z) 



(z) 



clM 



1 + 



0~M 



6 clM 



cm 



u 



- 1 



(4) 



In Eq. (O 6 c (z) = A c /D(z), where D(z) is the linear growth factor, 
cr M is the rms of primordial density fluctuations on the scale cor- 
responding to mass M, while S 3 (M) = fti^p 3 (M)lcr 4 u ls me nor " 
malized skewness. The third-order moment p 3 (M) can be written 

as 



H 3 (M) 



f Mr 



(*i)AW& 2 )M*(fc3)fla(ki,k2,k 3 )- 



dk 1 dk 2 d]i 3 

(2/r) 9 



(5) 



The function Mr(Jc) relates the Fourier transform of density fluctu- 
ations smoothed on some scale R to the relative peculiar potential, 
and is defined as 



2 T(k)k 2 
M R (k) = - 



3 H Q. m _o 



W R (k), 



(6) 



where T(k) is, in ou r analysis, the matter transfer function 
jEisenstein & HvJI 19981) and W R (k) is the top-hat window function. 
Eq|4] gives the correction to the PS mass function. However, we 
know that the PS mass function provides only an approximate fit to 
the results of N-body simulations. The commonly adopted proce- 
dure is then to assume that the same correction appearing in Eq.© 
can be applied to the best-fit Gaussian mass function to derive an 
accurate expression for the non-Gaussian one: 



n(M,z) = n (G) (M,z) 



n NG (M,z) 
n PS (M,z) 



(7) 



Here n { Q } {M, z) is the mass func tion in the refer e nce G aussian 
model computed according to the lSheth & Tormenl J2002T) recipe, 
while n NG (M,z) and n K {M,z) represent the IPress & Schechtej 
(1974) mass functions in the non-Gaussian (Eq. |4) and reference 
Gaussian model respectively. 

As already discussed, different prescriptions can be found in 
the literature that give different expre ssions for the mass fu nction 
in non-Gaussian models. For instance. Matarrese et alj (2000) used 
the saddle point approximation to compute the probability distri- 
bution of threshold crossing, and then truncated the resulting ex- 
pression to the skewness. The resulti ng mass function is h owever 
in agreement with the one obtained bv lLoVerde et alj d2008l) . as are 
the other recipes found in the literature. 

iGrossi etai] 12009) have shown that these analytic expressions 
are in agreement with N-body cosmological simulations, provided 
the linear threshold for collapse is corrected for ellipsoidal density 
perturbations, according to A c — > A c y[q, with q = 0.75. We adopted 
this correction in our calculations. 
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2.2 Bias 

The halo bias acquires an extra scale de pendence due to primordia l 
non-Gaussianity, that can be written as dMatarrese & V erde 2008) 

b(M, z, k) = b {0 \M, z) + Ab(M, z, k), (8) 

where 

Ab(M, z, k) = [b iG> (M, z) - l] 6 c (z)r R (k) . (9) 

The term T R (k) encapsulates the dependence on the scale and on 
the mass, and can be written for a local bispectrum shape as 

^_ r dt?M R (.oPo(ox 

i R (k)cr R Jo 



r*(*) 



8n 2 M R 



dp M R ( V«) 



Pdk) 



+ 2 



(10) 



where a = £ 2 + k 2 + 2p£k. The linear bias in the reference Gaussian 
model can be written as 



b {G) (M,z) 



1 + 

2p 



aA c 



1 



(D(z)cr M ) 



2p 



\ 2 P 



(ID 



(D(z)o- M ) 2 " + (^A c )- 

according t o the prescriptions of IMo & White! dl996h : 
ISheth & TormerJ dl999h : IShefh et al.l d200ll) . and with parameters 
set to p = 0.3 and a = 0.75. 

In order to obtain agreement with the results of numerical sim- 
ulations it is necessary to correct the linear overden sity for col- 
lapse , this time according to A c — > A c q with q = 0.75 dGrossi et al.l 
2009). Semi-analytic results, with t his correction, are also in agree- 
ment with the numerical results o f |Pillepichetai1d2010h (see also 
IPesiacqueT et al. 2009). We adopted this correction in the remain- 
der of our calculations. For figures showing the scale, mass and 
redshift dependenc e of t his correction to th e linear bias we refer to 
iTaruva et alj J2008I) and lFedeli et"al] d2009h . 



3 FISHER MATRIX ANALYSIS 

The Fisher Matrix formalism can be used to understand how accu- 
rately we can estimate the values of a vector of parameters p for a 
given model from one or more data sets, under the assumpti on that 
all parameters follow a Gaussian distribution (e.g.. lCashll 19791) . 
The information Fisher Matrix (FM hereafter) is defined as 



F, 



a)8 : 



(12) 



\dp a dp P l 

where £, is the likelihood of an observable (in our case the number 
of galaxy clusters in a given redshift and mass range or the averaged 
power spectrum of the cluster distribution). 

3.1 Number counts 



Following the a pproac h of iHolder et alj |2001) and 
Maiu mdar & Mohn 12003), the Fisher matrix for the number 
of clusters, A'/.,,,, within the /-th redshift bin and m-th bin in 
observed mass M" b , can be written as 



1 or/3 



z 



dN Lm dN u „ 1 
dp a dpfi Nij, 



(13) 



where the sums over / and m run over redshift and mass intervals, 
respectively. With this notation, it is Mf* = M,h r (z), where M t h r (z) 
is defined as the threshold value of the observed mass for a cluster 
to be included in the survey. Due to the selection function of any 
X-ray flux-limited survey, the value of Mtta(z) depends on redshift 
(see Section 147X1 . Therefore, the number of mass bins and their 
extent in our analysis will change with redshift accordingly. 

We write the number of clusters expected in a survey having a 
sky coverage AQ with observed mass between Mf b and M? b , , and 
redshift between zi and zi+i as 



AO. 



i; 



f ZM dz dV ["''' 
J ZI dzdCl J M cb 

dMn(M,z)p(M oh \\M). 



dM°' 



(14) 



In the above equation dV/dz is the cosmology-dependent comov- 
ing volume element per unity redshift interval and solid angle, 
n(M, z) the mass function, i.e. the number density of clusters with 
true mass M at redshift z (see Section ETI l. 

As alr eady mentioned, we assu me in the following the ex- 
pression by ISheth & Tormenl dl999h for the Gaussian halo mass 
function. We remind here that other calibrations of the halo mass 
functi on from simulations have been presente d by several authors 
(e.g.. IJenkinsetal .11200 ll ; IWarren et al] 120061 : iTinker et al.l 120081: 
ICrocce et alj|2010l) . While using the best-calibrated mass function 
is in fact impor tant when deri ving cosmological constraints from 
real data (e.g.. IWu et alJuOlOf) , it has only a minor impact when 
deriving forecasts on cosmological constraints. Indeed, what mat- 
ters for the latter is the total number of clusters expected in a given 
cosmological model, which is far more sensitive to the choice of 
reference cosmological and nuisance parameters than to the details 
of the mass fun ction fitting functio n. 

Following [Lima &Hu (2005) we assign to each cluster with 
true mass M a probability p{M \\M) of having an observed mass 
M ob , as inferred from a given mass proxy. Under the assumption of 
a log-normal distribution for the intrinsic scatter in the relation be- 
tween true and observed mass, with variance cr 2 nM , the expression 
for the probability is 

nh expt-.^M *)] 



V( 27r °"toM) 



(15) 



where 



x(M ob ) 



InM" 



B M - InM 



(16) 



Here we allow the relation between true and observed mass to be 
characterized not only by an intrinsic scatter, but also by a system- 
atic bias in the mass estimate, whose fractional value is given by 
B M . By inserting Eq.d 15t into Eq.d 14t for the cluster counts in a 
given mass and redshift interval, we obtain 



dz-r-rp- dMn(M,z) 
J dzdLl J Q 



An 

~2 J„ ~~ dzdfl . 
x [erfc(x m ) - erfc(x m+1 )] 



(17) 



with x„, = x{M ob m ) and erfc(x) the complementary error function. 
We note that the possibility of factorising the sky-coverage outside 
the integration relies on the assumption that clusters in our survey 
are detected over the same area of the sky down to the survey com- 
pleteness limit. The above expression can be easily generalized to 
include the possibility of a "flux"-dependent sky coverage. 
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We remark that we neglect clustering contribution to the noise 
(i.e. cosmic variance). Indeed, the Wide and the Medium surveys 
(see Table[T} covers large enough area over which cosmic variance 
is negligible. This may not be case for the Deep survey, which how- 
ever, as we shall see, does not bring much cosmological constraints 
(see Section |4~Tt . 



3.2 Power spectrum 

In order to include in our analysis the information from 
the clustering of galaxy clusters, we follow the approach by 
Majumdar & Mohr (2004). We define the Fisher Matrix for the 
power spectrum of galaxy clusters as 



a/3 ■ 



1 

(2*) 



d In P d (k m , Zi) d In P d (k„„ Zi) 



vfJtiAk, (18) 



where the sums in / and m run over redshift and wave number 
k bins, respectively dTegmarklll997l : iFeldman et al .11 19941) . In the 
above equations P cl is the average cluster power spectrum calcu- 
lated within the given redshift interval, 



7 .) = h 2 



(19) 



This amounts to weight the cluster power spectrum, P° l (k, z), ac- 
cording to the square of the number density of clusters, N(z), 
that are included in the survey at redshift z. In turn, the cluster 
power spectrum P d (k, z) is expressed in terms of power spectrum, 
P(k, z), of the cosmic density fluctuations according to P cl (k, z) = 
t%ff(k, z) P(k, z); here the effective bias is defined as the linear bias, 
introduced in Sect. 2, weighted by the mass function, 



b eff (z,k) 



dMn(M, z) erfc[x(M thr )] b(M, z, k) 
£°dMn(M,z)eifc[x(M tbl )] 



(20) 



The bias parameter b(M,z,k) acquires the dependence on the 
wavenumber k predicted by non-Gaussian models (see Sect. 12.2b . 
Finally, the quantity V eif (k, z) in Eq.dl8t is the effective volume ac- 
cessible by the survey at redshift z at wavenumber k. This effective 
volume is weighted by the shot noise level l/7V(z), so that 



V eff (k,z) = V (z) 



N(z)P c ,(k,z) 



\+N(z)P c £k,z) 



(21) 



with V a {z) the total comoving volume covered by the redshift bin 
centred on z- In this way, constraints at redshift z are mostly con- 
tributed by wavemodes k, which maximize N(z)P c i(k, z) and make 
V e ff approach V . 

An alternative approach to include clustering in formation in 
derivin g FM survey foreca sts has been proposed by iLima & 53 
d2005h and applied also bv lQguri| [2009) for constraints on non- 
Gaussian models from cluster surveys. In this approach, one makes 
a partition of the sky area covered by a survey into regular cells of a 
fixed angular size and then computes the fluctuations in the cluster 
counts within such cells. Since this method does not explicitly in- 
clude the covariance between counts within different cells, it only 
samples clustering at a fixed angular scale (i.e. at a single physical 
scale for a fixed redshift). On the other hand, extracting the full in- 
formation content in the scale dependence of the power spectrum 
is quite important when constraining non-Gaussian mo dels, whose 
unique signature is given by the scale-dependent bias. lCunha et al.l 
(2010) used the count-in-cell approach by also including the in- 
formation from the covariance. Therefore, the information on the 



large-scale power spectrum, in the count-in-cell approach, is con- 
veyed by the covariance terms. In our approach the different scales 
are weighted by the effective volume, defined by Eq.d21b. 

In our analysis we assume the following reference values for 
the cosmolog ical parameters, con sistent with the WMAP-7 best- 
fitting model dKomatsu et al.|[20~ich : fl„, = 0.28 for the present-day 
matter density parameter, <x 8 = 0.81 for the normalization of the 
power spectrum, flj. = for the contribution from the curvature, 
w(a) = wo + (1 - a) w a with wo = -0.99 and w a = for the 
Dark Energy equation of state, Q. h = 0.046 for the contribution of 
baryons to the density parameter, h = 0.70 for the Hubble parame- 
ter, n = 0.96 for the primordial spectral index and /nl = for the 
non-Gaussianity parameter. Therefore, we have in total 9 cosmo- 
logical parameters, which are left free to vary in the computation of 
the number counts and power spectrum Fisher Matrices of Eqs.dl3t 
andl[T8j. 

In the following, all the results presented are based on 
adding the Fisher Matrix for the Planck CMB experiment to 
those from the cluster surveys. We derive the cosmological con- 
straints from P l anck f ollowing the description laid out by the DETF 
lAlbrecht et al] d2009l) and use the method described in Ras sat et"al] 
(2008). We conservatively assume that we will only use the 143 
GHz channel as science channel. This channel has a beam of 
#fwhm =7.1' and sensitivities of a T = 2.2/uK/K and o> = 4.2p.K/K. 
We take / sky = 0.80 as the sky fraction in order to account for galac- 
tic foregrounds. We use as a minimum f-mode, £ min = 30 in order to 
avoid proble ms with polarization foregrounds. As described in the 
DETF report dAlbrecht et alj2009l) we choose as fiducial parameter 
set 6 = (a) m , 9 S , In A s , u> b , n s , t), where S is the angular size of the 
sound horizon at last scattering, In As is the logarithm of the pri- 
mordial amplitude of scalar perturbations and t is the optical depth 
due to reionization. After marginalization over the optical depth 
we then calculate the Planck CMB Fisher matrix in the parame- 
ters (£!„,, O rfc , h, cr 8 , £i b , w , w a ,n s ) by using th e appropriate Jaco - 
bian of the involved parameter transformation dRassat et al]|2008l) . 
We point out that the Planck FM is computed for Gaussian per- 
turbations. Therefore, while it adds quite strong constraints on the 
other cosmological parameters, especially on the curvature, it does 
not add any constraints on /nl. 



4 ANALYSIS AND RESULTS 
4.1 Characteristics of the surveys 

The commonly adopted procedure to estimate cosmological fore- 
casts for X-ray surveys is based on calibrating a flux limit for clus- 
ter detection (which generally corresponds to the detection of few 
tens of net photon counts in an extended source), and to use fluxes 
as proxies to cluster masses. In order to account for the uncertain 
knowledge of the relation between X-ray luminosity and mass, sev- 
eral authors h ave proposed to follow the so-cal l ed self-calibration 
metho d (e.g.. iMaiumdar & Mohri [20031 12004| ; ILima & Hij |2004[ 
2005). In this approach the relation between mass and observable 
is defined up to an intrinsic scatter, also including the possibility of 
a systematic bias in the estimate of cluster masses. The parameters 
defining the relation between mass and observable are then treated 
as fitting "nuisance" parameters to be determined along with the 
relevant cosmological parameters. 

An alternative approach, adopted to derive cosmological con- 
straints from cluster surveys, is, instead, based on restricting the 
analysis to a relatively small subset of galaxy clusters for which 
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Table 1. Characteristics of the X-ray surveys (see also lGiacconi et al . 2009. 
Column 2: sky coverage fl (in sq.deg.); Column 3: flux limits for detec- 
tion of extended sources in the [0.5-2] keV energy band (units of 10~ 14 erg 
s _1 cm~ 2 ); Column 4: flux limits defining the bright subsamples (see text). 





n 




fbr 


Wide 


20000 


0.5 


15.0 


Medium 


3000 


0.1 


3.0 


Deep 


100 


0.01 


0.3 



deeper X-ray follow-up observations provide measures of mass 
proxies which are more closely relate d to cluster masses (e.g., 
IVikhlinin et"al]|2009l : iMantz et alj|2009l) . Examples of such robust 
mass proxies are the so-called Y x = M gas T x , define d as the product 
of ga s mass and temperature at a given radius (e.g.. lKravtsov et ahl 
2006), or M gas itself. They are robust in the sense that, based on 
simulations, their relation with cluster mass has an intrinsic scatter 
of 10 per cent or less, especially when emission from cluster cores 
can be resolved and/or removed. While having obvious advantages, 
the limitation of relying on such low-scatter mass proxies is that 
they can be measured for a relatively small number of clusters. In 
this sense, one has to compromise between the requirement to keep 
under control the systematics in mass measurements and the need 
to cover a relatively large redshift baseline with adequate statistics. 

Ideally, one should carry out an X-ray survey at such a good 
sensitivity that low-scatter mass proxies can be measured in survey 
mode for a large number of galaxy clusters. This would allow an 
"educated" self-calibration analysis of the cluster surveys, in which 
the systematics in mass measurements can be kept under control. 

In the following, we derive forecasts for three X-ray surveys, 
which are complementary in terms of sensitivity and sky coverage, 
inspired to the survey strategy devised for the Wide Field X-ray 
Telescope (WFXlQ), recently proposed to the Astro-2010 Decadal 
Survey panel (e.g.. iMurray & WFXT Teamll2010l : iGiacconi et al.l 
2009; Vikh linin et alj|2009h . This telescope combines a large col- 
lecting area and field of view with a sharp PSF over the entire field 
of view. Thanks to these characteristics, the WFXT has the poten- 
tial of detecting a large number of clusters out to z ~ 2 and of 
measuring Y x or M sas with good precision down to relatively low 
fluxes. While describing the characteristics of the WFXT is outside 
the scope of this paper, we point out that such an instrument would 
have the capability of measuring low-scatter mass proxies down to 
cluster fluxes which are comparable to fluxes at which clusters are 
just detected in curre nt and future X-ray telescopes. Based on the 
results presented by IGiacconi et aD §009), we give in Table [T] the 
limiting fluxes (in the 0.5-2 keV energy band) at which WFXT will 
detect a cluster as an extended source in three surveys: a Wide Sur- 
vey, which will cover all the extragalactic sky (20000 deg 2 ) with 
a sen sitivity ~ 500 tim es better than the ROSAT All Sky Survey 
(e.g.. lVogesetIfl [l999h: a Medium Survey, which will reach over 
3000 deg 2 flux limits compara ble to those of the deep Chandra and 
XMM deep COSMOS fields dCappelluti et~ai1l2009n : a Deep Sur- 
vey that will reach over 100 deg 2 a sensitivity similar to those of the 
deepest Chandra pointings. Thanks to the large collecting area and 
field of view of WFXT, and taking advantage of its good angular 
resolution (5 arcsec half energy width, approximately constant over 
the whole field-of-view), the se surveys could be completed within 
a five-year mission duration dMurrav & WFXT Teamll201o|) . 



In principle, a unique flux limit is not sufficient to define a 
completeness criterion in an X-ray survey. In fact, due to vignetting 
and PSF variation with off-axis angle, the flux limit for the detec- 
tion of a source at a given signal-to-noise varies across the field of 
view. For this reason, rather than a flux limit, one should calibrate 
a flux-dependent sky coverage. Owing the approximate uniformity 
of the WFXT PSF, we expect such a sky coverage to be quite steep 
around the flux limits reported in TableQ] so that we ignore its flux- 
dependence in the following analysis. In order to convert these flux 
limits into mass limits, we u se the relation be tween X-ray lumi- 
nosity and M500 calibrated bv lMaughanl J2007t) . where masses are 
recovered from Y x , using Chandra data for 115 clusters in the red- 
shift range 0. 1 < z < 1.3. Among the fitting expressions reported in 
Table 1 of that paper, we choose the relation between L x and M500, 
obtained without excising the core region within 0. 15 Rsm' 



CE(z) 



M 5 , 



4 x 10 14 M e 



(22) 



withC = 5.6, B = 1.96. 

The reason for this choice is that we did not attempt to model 
the core contribution in computing the flux limits reported in Table 

m 

We show in Figure Q] the redshift dependence of the limit- 
ing mas.<0 Mjoo associated to the survey flux limits. The value of 
the virial mass, which is the relevant quantity entering in the mass 
function choosed in ou r analysis, is obtained from M500 , following 
iHu & Kravtsovl 1 2003). by adopting the NFW halo density profile 
dNavarro et al J 19971) for the reference cosmological m odel, assum- 
ing c = 5 for the concentration parameter (see also Shang et al. 
2009). Also shown in Figure Q] with the green short-dashed curve 
is the mass limit corresponding to a flux 30 times brighter than the 
flux limit for cluster detection in the Deep Survey. Since cluster 
detection in the Deep survey corresponds to about 200 net counts, 
such a mass limit is for clusters for which about 6000 counts would 
be available. With such a large number of counts one has a precise 
measurement of robust mass proxies such as Y x and M gas . We point 
out that clusters identified in the Deep Survey allows one to cali- 
brate mass proxies down to fluxes which are lower than the flux for 
cluster identification in the Wide Survey. Such brighter flux lim- 
its define the bright samples of clusters, for which a direct mea- 
surement of a mass proxy can be carried o ut within the sa me sur- 
veys. Extrapolating the M-L x relation bv lMaugharj J2007h at faint 
fluxes of our surveys would imply unrealistically small mass limits 
at low redshift. For this reason, we decided to use a lower limit of 
M500 = 5 x 10 13 /r'M in the definition of the selection function 
(shown with the horizontal dotted line in Figure [TJ. In fact, this is 
comparable to the lowest m ass down to which m ass proxies have 
been calibrated so far (e.g.. IVikhlinin et af. 2009). We also point 
out that, for the sake of simplicity, we do not include the cosmol- 
ogy dependence of the selection function in our analysis. 



4.2 Nuisance parameters 

Besides the nine cosmological parameters, our FM analysis should 
also constrain the nuisance parameters which specify the redshift 
dependence of the fraction mass bias B M and the intrinsic scatter 
finM (in our analysis we do not consider the case of a possible 



1 http://www.eso.org/~prosati/WFXT/Overview.html, 
http://wfxt.pha.jhu.edu/ 



2 Here and in the following we indicate with Ma the mass contained within 
the radius encompassing an average density of A times the critical cos- 
mic density p cr . 
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Figure 1. The redshift dependence of the M500 mass thresholds, corre- 
sponding to the flux limits for cluster detection for the three surveys, as 
reported in Table 1 : the Wide, Medium and Deep surveys are shown with 
the solid (red), the dotted (blue) and dot-dashed (cyan) curves, respectively. 
The short-dashed (green) curve corresponds to a flux-limit which is 30 times 
brighter than the flux limit for cluster detection in the Deep survey. The hor- 
izontal dotted line marks the mass limit of M500 = 5 X 10 13 /i~'M Q below 
which we discard clusters in our analysis. 

mass dependence of these parameters). According to iLima & Htj 
(2005), we assume the following parametrization for such redshift 
dependencies: 

B M (z) = B MSS {\+z) a 
o"i„m(z) = ar lnM0 (l + zf ■ (23) 

In this way, we have four nuisance parameters, B M .o, cinM.o, ct 
and B. A negative value for B M corresponds to a mass under- 
stimate and, therefore, to a smaller number of clusters included 
in a survey, for a fixed selection function. The presence of the 
mass bias accounts for the possibility of a violation of hydro- 
static equilibrium in the estimate of X-ray masses, on which the 
observable-mass scaling relation is calibrated. A number of inde- 
pendent analyses of a variety of cosmological hydrodynamic sim- 
ulations of galaxy clusters converge to indicate that hydrostatic 
mass estimators provide unde restimates of true mass within R 500 
by about 10-15 per c ent (e.g..lRasia et alJl2005l:lNagai et alj|2007l : 
lAmeglio et ai]|2009l ; IPiffaretti & Valdarninill2008l) . Quite reassur- 
ingly, such results also agree with the observational results on the 
comparison between clust er masses estimated with weak lensing 
and with X-ray data (e.g., Mah davi et al J 1 20081) . In the following 
we assume B M0 = -0.15 as a reference value for the mass bias 
and regarding its evolution, we take a = as a reference value. As 
for the intrinsic scatter, it has the effect of increasing the number 
of clusters included in the survey. In fact, the number of low-mass 
clusters that are up-scattered above the survey mass limit is always 
larger than the number of rarer high-mass cl usters which are down- 
scattered below the same mass limit (e.g.,[Cunha 2009, and ref- 
erences therein). As a reference value, we assume cr\ nM ,o = 0.25, 
consiste nt with the instrin sic scatter in the M 500 -L x relation mea- 
sured bv lMaughan (2007), with B = for its evolution. We stress 
here that, following iTima & Hul J2005h . we use the variance tr. 



and not the scatter as the parameter to be varied in our Fisher matrix 
analisys. In fact, this quantities controls the excess of up-scattered 
and down-scattered clusters with respect the total number. 

In the following we will not assume any prior for these four 
nuisance parameters in our reference analysis. We will refer to it as 
the no prior analysis. 

On the other hand, already available data allow one to set 
constraints on the value of the mass bias from the comparison be- 
tween X-ray and lens ing cluster mass measurements. For instance, 
Vikhl inin et all d2009h compared weak lensing and Chandra X-ray 
mass measurements for a rather small sample of low-z clusters and 
concluded that the mass scale can already be calibrated with a sta- 
tistical unce rtainty of about 10 p er cent. A similar result has been 
obtained by IZhang et alj J2010f) from the comparison of XMM- 
Newton X-ray masses and weak lensing masses for a set of 12 
nearby clusters. Owing to these results, we also consider the case of 
the mass bias parameter to be known with a precision AB M0 = 0.05. 
An improvement by only a factor of two with respect to the present 
in the calibration of the cluster mass scale is probably overconser- 
vative, owing to the orders-of-magnitude increase in the number of 
clusters with precise X-ray mass measurements to be provided by 
the three surveys and precise lensing mass measurements from both 
ground-based and space telescopes. As for the evolution, we as- 
sume Ao- = 1 as a prior, which would correspond to an uncertainty 
in the mass bias calibration at 1 = 1 comparable to that calibrated 
at present for nearby clusters. Regarding the prior on the intrinsic 
scatter, we assume A<xi„ M ,o =0.1 and AB = 1 . We expect these to be 
rather conservative choices, in view of the large number of clusters 
that should be made available by future X-ray and optical/near-IR 
surveys at both low and high redshift. In the following, we refer to 
these choices of the priors for the nuisance parameters as the weak 
prior analysis. 

Finally, we also refer to the strong prior analysis for the case 
in which nuisance parameters are assumed to be known with little 
uncertainty and are thus kept fixed to their reference values. While 
this assumption is expected to be unrealistic for fluxes reaching the 
limiting flux for cluster detection in the three surveys, it may be 
rather plausible in case we restrict the analysis to the bright sub- 
samples. 

While we use the no prior choice as a reference for the 
analysis of the surveys down to the detection limit, we will show 
how much we would gain in terms of constraining power by 
using instead the weak prior and the strong prior assumptions. 
Finally, we will present results for the bright surveys for the case 
of strong priors on the nuisance parameters. This will allow us to 
judge the trade-off between having statistically richer (i.e. lower 
flux limit) surveys with a less controlled calibration of the cluster 
mass-luminosity conversion, and a brighter flux limit with better 
controlled nuisance parameters. 

As a first check, we verify that our reference model provides 
a good fit to present data. We show in Figure[2]the comparison be- 
tween the cluster flux number counts observed in the [0.5-2] keV 
energy band from the ROSAT Deep Cluster Survey dRosati et al.1 
1998,'l2002j) and the prediction of our reference model. This is not 
surprising, owing to the fact that the reference cosmological pa- 
rameters agree with co nstraints based on the evolution of the clus- 
ter mass function (e.g., iBorganiet al]|200ll : I Vikhlinin etaTi r2009: 
Mantz et al. 2009). The good agreement between available obser- 
vational data and predictions of our reference model indicates that 
the latter can be used to provide a realistic extrapolation of the evo- 
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Figure 2. The comparison between observed cumulative cluster flux num- 
ber counts (symbols with errorbars) and predictions from the reference 
model (dotted curve, see text). Observational results refer to counts in 
the [0 .5-2] keV band from the ROSAT Deep Cluster Survey teosati et alj 
|2002|) . with errorbars corresponding to lcr Poissonian uncertainties. 

lution of the cluster mass function over redshift and mass ranges 
which are not probed by currently available data. 

We show in Figure [3] the cumulative redshift distributions for 
the clusters to be detected in the three surveys (left panel) and for 
the bright samples (right panel). Overall, the three WFXT surveys 
would yield about 3x 10 6 detected clusters, out of which ~ 7.5 x 10 4 
clusters should be found at z > 1. This will provide an improve- 
ment by about four orders or magnitude with respect to the ~ 10 
z > 1 clusters currently confirmed. At the same time, we expect 
to have about 2 x 10 4 clusters with robustly measured mass prox- 
ies, of which ~ 4000 would lie at z > 0.5. This would increase 
by more than two orders of magnitude the number of clusters for 
which mass proxies have been measured above this redshift, af- 
ter intensive follow-up Chandra observations of clusters identified 
in ROSAT-based surve ys (e.g., IVikhlinin et al.l l2009; Man tz et al.l 
l2009tlEttori et al.l2009l) . We stress that, despite the small area cov- 
ered, the Deep survey provides the dominant contribution to the 
bright sample at z > 1. This highlights the important role that the 
Deep survey has in providing mass proxies at high redshift. Al- 
though predicting the number of extremely distant clusters is highly 
uncertain, owing to the unknown evolution of the mass luminosity 
relation above z — 1, we foresee that ~ 10 3 clusters would be de- 
tected at z~ 2, with mass measurements available for few tens of 
them. In order to quantify the increase in sensitivity provided by the 
WFXT surveys with respect to currently planned X-ray missions, 
we computed the redshift distributions expected for the surveys to 
be carried out by the eROSITA satellit^. In the three years of oper- 
ation eROSITA is expected to find, for our reference model, = 4000 
clusters at z > 1, while virtually no clusters with precise measure- 
ments of mass proxies would be found at z > 0.5 in the bright 
surveys having flux limits 30 times higher than for detection. 

An important lesson learned from the Chandra follow-up ob- 

3 http://www.mpe.mpg.de/erosita/MDD-6.pdf 



servations of distant clusters is the fundamental relevance of a sharp 
PSF to excise the contribution of cool cores in the measurements 
of mass proxies . Indeed, as show n by a number of authors (e.g., 
Maughai] |2007l; IPratt et alj[201oh . excising the core contribution 
suppresses by a substantial amount the intrinsic scatter in the rela- 
tion between cluster X-ray observables and masses. The sharp PSF 
with which the WFXT would carry out the surveys, should guaran- 
tee a good control of the cool core contribution in the calibration 
of mass proxies, even to z ~ 1.5 without the need of follow-up ob- 
servations with other higher resolution X-ray telescopes. We also 
emphasize that a sharp PSF, such that expected for WFXT, has the 
additional benefit of easing the subtraction of the point source con- 
tribution when measuring cluster fluxes. This allows an accurate 
assessment of the survey flux completeness. 

4.3 Constraints on non-Gaussianity 

Having defined the reference cosmological model and mass-flux 
conversion, we present now forecasts on constraints from non- 
Gaussian models. These results will be shown in terms of con- 
straints on the cr 8 -/ NL plane after marginalizing over the other cos- 
mological and nuisance parameters. The reason for this choice is 
that, for a fixed Friedmann background, <r 8 and / NL are the two 
parameters which determine the timing of structure formation and, 
therefore, the evolution of number density and large-scale cluster- 
ing of galaxy clusters. In the following, we will always show con- 
straints on the o"8-/nl plane at the 68 per cent confidence level. 

In analogy with the figure-of-merit introduced to quantify the 
constraining pow er of an experiment for the Dark Energ y equa- 
tion of state (e.g.. lAlbrecht et al.ll2006tl2009l;IWang| l2008), we in- 
troduce a figure-of-merit for the timing of structure formation in 
non-Gaussian models: 

FoMsft = (det[Cov(o- 8 ,/ NL )])- 1/2 , (24) 

where Cov(os, /nl) is the covariance matrix between <t$ and /nl, 
which is obtained by inverting the FM and marginalizing over all 
the other parameters. 

In the computation of the FM of Eq.lll3t. the /V;,„ number 
counts are computed out to z = 2 within 60 constant redshift bins 
and within observed mass bins having width A log M - 0.1, extend- 
ing from the lowest mass limit determined by the selection function 
(see Fig.QJ and an upper mass limit of 10 16 /j~'M o . We verified that 
a finer binning does not add more information in the Fisher Matrix. 
As for the redshift bins, their size is larger than the redshift errors 
obtainable from optical spectrscopy, and comparabl e to the typical 
uncer tainties from red-sequence redshift estimates (Gladd ers et al.l 
|2007|) . 

As for the computation of the FM for the power spectrum 
of Eq.fTgJ, we used wavenumbers in the range k max > k > 
0.001 Mpc~', independent of redshift. Using an arbitrary small 
value of k does not change the final results, since extremely large 
wave modes are not sampled by the surveys and, therefore, do not 
provide any contribution to the FM. As for the value of k max , the 
chosen value should represent a compromise between the needs of 
maximising the amount of information to be extracted from the sur- 
vey and of avoiding the contribution of small-scale modes where 
the validity o f the linear bias model is compromised by the onset 
non-linearity dPercival & Whitd2009l) . In the following, we will as- 
sume k max = 0.3 Mpc~', and will also show the sensitivity of the 
results to this choice. 

As for the redshift binning, the average cluster power spec- 
trum defined by Eq .d 19t is computed by integrating over redshift 
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Figure 3. The cumulative redshift distribution for the three surveys. The left panel is for all clusters to be detected down to the survey flux limits, while the 
right panel is for the clusters in the bright samples, corresponding to a 30 times higher flux limit. In both panels solid (red), dotted (blue) and dot-dashed (cyan) 
curves represet the Wide, Medium and Deep surveys, respectively, while the short-dashed (green) curve represents the sum of the three. 



intervals having constant width Az = 0.2. This coarser binning, 
with respect to that used for the analysis of number counts, is dic- 
tated by a compromise between the need of extracting the maxi- 
mum amount of information from the clustering evolution and the 
request of negligible covariance between adjacent z-intervals (e.g., 
IStriletaLlk oiO). Indeed, the contribution from different z-bins can 
be added in the defining of FM of Eq.J18t only if they carry statis- 
tically independent information. 

We show in Figure [4] the constraints on the /nl and erg pa- 
rameters computed from the number counts and from the power 
spectrum within the Wide survey, by assuming no prior on the val- 
ues of the nuisance parameters. This plot clearly demonstrates the 
strong complementarity that number counts and large-scale cluster- 
ing have to constrain erg and /nl: while number counts are highly 
sensitive to the value of <x 8 , the weak sensitivi ty of the high-end 
of the mass function to non-Gaussianity (e.g.. iFedeli et alj|2009l 
and references therein) provides only very weak constraints on /nl; 
conversely, the scale-dependence of bias makes the power spectrum 
a powerful diagnostic for non-Gaussianity, while providing only 
loose constraints on erg. 

If we combine all the information obtainable from the three 
surveys, we obtain the constraints shown in Figure [5] Most of the 
constraining power is provided by the Wide survey, with only little 
information on structure formation timing carried by the Medium 
and Deep surveys. There are two main reasons for this. First, the 
Wide survey provides the largest statistical baseline out to z — 1.5, 
when including all clusters down to the mass limit corresponding 
to detection (see left panel of Figure [3]l. This implies a better de- 
termined mass function and, therefore, stronger constraints on erg. 
Second, the larger area coverage of the Wide survey allows it to bet- 
ter sample long-wavelength modes, where the scale-dependence of 
the bias induced by non-Gaussianity can be better assessed, thus 
turning into stronger constraints on / NL . As shown in Table [2] the 
value of FoM SFT for the combination of the three surveys is in fact 
dominated by the Wide Survey. 

As already mentioned, we have assumed in our analysis k max = 
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Figure 4. Constraints at the 68 per cent confidence level on non-Gaussian 
parameter /nl and power spectrum normalization erg coming from number 
counts alone (short-dashed green curve), power spectrum alone (dotted blue 
curve) and from the combination of the two (solid red curve). The analysis 
refer to the clusters detected in the Wide Survey. We marginalized over all 
the other parameters. No prior is assumed for the values of the nuisance 
parameters. The Fisher Matrix from Planck experiment is included in the 
calculation of all constraints. 

0.3 Mpc~'. In order to quantify the sensitivity of our results to the 
adopted k max value, we show in Figure|6]how the constraints change 
for k max = 0.1 Mpc~' and k max = IMpc -1 . The smaller value is 
close to the scale of non-linearity at z = 0, although it is prob- 
ably too conservative at high redshift, z ~ 1. As expected, de- 
creasing k max makes the constraints slightly looser, due to the lower 
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Table 2. Figure-of-merit of structure formation timing, FoMsft [see Eq. 
1241 . and r.m.s. uncertainty in the non-Gaussian parameter, (Tf NL , for the 
three surveys, and for their combination, assuming different priors for the 
nuisance parameters. Columns 3-6 show the results for the Wide, Medium 
and Deep Surveys, and for the combination of the three. 







Wide 


Medium 


Deep 


Total 


Detection - no pr. 


FoMsft 


33.1 


8.5 


0.4 


39.2 




°7nl 


11.3 


18.5 


84.2 


10.4 


Detection - weak pr. 


FoMsft 


33.3 


8.8 


0.6 


39.4 




°7nl 


11.3 


18.4 


84.3 


10.4 


Detection - strong pr. 


FoMsft 


157.2 


49.3 


3.0 


183.2 




°7nl 


11.2 


18.0 


80.9 


10.3 


Bright - strong pr. 


FoMsft 


7.3 


14.3 


3.0 


22.3 




°7nl 


55.9 


45.9 


85.7 


33.8 
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Figure 5. Constraints at the 68 per cent confidence level on non-Gaussian 
parameter /nl and power spectrum normalization erg from the Deep, 
Medium and Wide surveys (dot-dashed cyan, dotted blue and solid red 
curves, respectively), by combining number counts and power spectrum in- 
formation, by using no priors on the nuisance parameters. Also shown with 
the short-dashed green curve are the constraints obtained from the combina- 
tion of the three surveys. No prior is assumed for the values of the nuisance 
parameters. The Fisher Matrix from Planck experiment is included in the 
calculation of all constraints. 



amount of information included in the Fisher Matrix of Eq.dlSt. 
Correspondingly, the value of the figure-of-merit decreases from 
FoMsft = 39.2 to 33.1, with the uncertainty on /nl increasing only 
from <t Al = 10.4 to 12.1. Increasing instead k mal to 1 Mpc~' does 
not lead to any significant improvement of the constraints. In fact, 
given the level of Poisson noise associated to the cluster distribu- 
tion, high frequency modes are not adequately sampled and, there- 
fore, adding them to the analysis does not add significant informa- 
tion. 

The contribution of information to the Fisher Matrix carried 
by the power spectrum at different redshifts and wavenumbers can 
be understood by looking at the dependence of the effective vol- 
ume, V e ff, on the power spectrum, which is set by the bias param- 




0.82 



Figure 6. Constraints at the 68 per cent confidence level on non-Gaussian 
parameter /„/ and power spectrum normalization o-g from the combination 
of the three surveys, when changing the maximum value of the wavenum- 
ber, k max , for which power spectrum information are included in the Fisher 
Matrix of Eq.{T8}. Solid (red), dotted (blue) and short-dashed (green) curves 
correspond to k max = 1, 0.3 and 0.1 Mpc~', respectively. No prior is as- 
sumed for the values of the nuisance parameters. The Fisher Matrix from 
Planck experiment is included in the calculation of all constraints. 



eter, and on the level of Poisson noise, which is set by the number 
density of clusters. Following Eq. d21l >. we define the quantity 



W v (k,z) 



N(z)Pci(k,z) 



\+N(z)P a (k,z) 



(25) 



which gives the weight carried by the wavenumber k to the compu- 
tation of the clustering Fisher Matrix at redshift z. In the left panel 
of Figure [7] we show the redshift dependence of the effective vol- 
ume computed within redshift intervals of constant width Az = 0.2, 
for different values of k, and compare them to the total comoving 
volume computed within the same redshift intervals. The effective 
volume lies always well below the total comoving volume: this is 
the consequence of the relatively low value of the cluster number 
density, which makes Poisson noise always dominating. While the 
total comoving volume V increases with redshift, the effective vol- 
ume V e ff starts declining after reaching a maximum, at z — 0.5, 
for all wavenumbers. As for the dependence on k, at a fixed red- 
shift, the value of V e ff decreases for both very high and very low 
values of k. As shown in the right panel of Figure UJ the value of 
the weight function W v (k, z) is maximized at k 0.01 Mpc~'. In 
fact, for W,, « 1 (i.e. N(z)P c i(k,z) «: 1), the ^-dependence of 
W v reflects that of P c /. Poisson noise is, again, responsible for the 
low values of W v , well below unity. Decreasing of the level of this 
noise would require increasing the number density of objects to be 
included in the survey. This could be accomplished in principle by 
decreasing the mass threshold. However, this would require bring- 
ing into the surveys low-mass clusters and groups, for which our 
parametrization of the mass-observable relation may not still be 
valid. 

So far, we presented results by assuming prior on cosmolog- 
ical parameters from Planck experiment and no prior knowledge 
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Figure 7. Left panel: the redshift dependence of the effective volume, defined as in Eq. j2U . within redshift intervals of constant width Az = 0.2, for four 
values of the wavenumber k. Short-dashed (green), solid (red), dotted (blue) and dot-dashed (cyan) curves correspond to k = 0.001, 0.05, 0.1 and 0.2Mpc~', 
respectively. Right panel: the dependence on the wavenumber of the weight W,,(k, z), defined as in Eq. j25t . at three different redshifts. Solid (red), dashed 
(green) and dotted (blue) curves are for z = 0.5, 1 and 1,5, respectively. 



on the nuisance parameters. We want to stress that, as already dis- 
cussed in Sect. 14.21 this is probably too much a conservative ap- 
proach, in view of the calibration of the relation between robust 
mass proxies (e.g., Y x and M gas ) and X-ray luminosity for a large 
number of clusters within the planned surveys. In Figure[8]we show 
the effect of assuming a prior knowledge of the nuisance parame- 
ters. If we assume the weak priors for these parameters (see Sect. 
14.2b . constraints are only slightly improved. Quite interestingly, 
even assuming the strong prior (i.e. nuisance parameters fixed) im- 
proves the constraints on cr 8 , while having a smaller impact on 
those for /ml- Indeed, we find that the error on non-Gaussianity 
only decrease from o-f NL a 10.4 to 10.3 when passing from the 
no prior to the strong prior assumption, while the figure-of-merit 
increases from FoMsft — 39.2 to = 183.2 (table:[2}. 

To better understand the reason for the weak dependence of 
the / NL constraints on the uncertain knowledge of the nuisance pa- 
rameters, we show in Figure [9] by how much number counts and 
effective bias change with respect to the value that they take in 
the Gaussian case, as we vary the mass bias parameter B M (left 
panels) and the intrinsic mass-scatter o"i nM (right panels). As a 
reference value for the non-Gaussianity, we take here / NL = 10, 
which is comparable to the forecasted precision with which non- 
Gaussianity can be constrained from our analysis. At z = 0.5 the 
deviation of the number counts from the non-Gaussianity (upper 
panels) varies only by about one part over thousand when a gener- 
ous range of variation is allowed for both B M and cr\ nM , with only 
a slighty higher sensitivity to these parameters at z = 1. In the 
bottom panels of Figure [9] we show the sensitivity of the effective 
bias on nuisance parameters for different values of the wavenum- 
ber k. Results are shown at z = 0.5 which is close to the redshift 
where the effective volume V e ff reaches its maximum value (see 
left panel of Figure [7}. For the level of non-Gaussianity assumed 
here, the deviation from the Gaussian effective bias is negligible at 
the wavenumbers, k a 0.01 Mpc -1 , which are mostly weighted in 
the computation of the Fisher Matrix (see right panel of Figure[7](. 
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Figure 8. Constraints at the 68 per cent confidence level on the non- 
Gaussian parameter /nl and power spectrum normalization <t% by assuming 
no prior (solid red curve), weak prior (dashed green) and strong prior (dot- 
ted blue) on the nuisance parameters. All constraints are obtained by com- 
bining cluster number counts and power spectrum information for the three 
surveys together. The Fisher Matrix from Planck experiment is included in 
the calculation of all constraints. 

However, as expected, the effect of non-Gaussianity on b e ff shows 
up at very large scales, with a deviation with respect to the Gaus- 
sian result by ~Z 40 per cent for k a 10~ 3 Mpc -1 . This highlights the 
importance for future surveys to have a highly uniform calibration 
of the selection function over large area of the sky, for them to be 
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Figure 9. Effect of changing the values of the nuisance parameters of mass 
bias Bu (left panels) and intrinsic mass scatter <T\ nM on the deviations of 
number counts and effective bias from the Gaussian case. The results shown 
here are obtained by assuming a mass-limit of 10 14 M o . The reference value 
of non-Gaussianity assumed here is /nl = 10. N is the ratio between the 
redshift distributions in the Gaussian and non-Gaussian cases, while B is the 
ratio between the effective bias, as defined in equation 1201 in the Gaussian 
and non-Gaussian cases. Upper panels show the results for number counts 
and effective bias at two different redshifts, z = 0.5 and z = 1, Lower panels 
show results for the effective bias at z = 0.5 for different wavenumbers, 
k = 0.1, 0.01 and 0.001 Mpc~'. 



able to appreciate any subtle scale dependence of the bias parame- 
ter. Also in this case, any variation with the value of the nuisance 
parameters is far smaller than the deviation from Gaussianity. This 
justifies the weak dependence of the / NL constraints on the uncer- 
tain knowledge of the cluster mass calibration. 

As already mentioned in Sect. 14. II the bright surveys contain 
only clusters for which a precise estimate of robust mass proxies, 
such as Yx and M gas , can be obtained. In this way, the analysis of 
the bright subsamples offers us the possibility of testing the trade- 
off between having smaller samples, for which strong prior on the 
nuisance parameters can be assumed, and of larger samples with a 
less certain knowledge of such parameters. In Figure [TO] we show 
the constraints expected for the bright subsamples of the three sur- 
veys. We note that this time the role of the Medium and of the Wide 
surveys are reversed, with the former providing the most stringent 
constraints. This is not surprising, since the number of bright clus- 
ters in the Wide Survey falls below that of the Medium Survey al- 
ready at z — 0.1 (see Figure [3). Therefore, the richer statistics of 
distant clusters in the Medium Survey more than compensates the 
more efficient sampling of small-A: modes offered by the Wide sur- 
vey. In terms of figure-of-merit, we note that the resulting value 
from the combination of the three bright surveys is FoM S ft — 22.3: 
this implies that assuming strong priors for nuisance parameters in 
the bright surveys has less constraining power than using no priors 
for the surveys defined down to the detection flux limit. 



5 DISCUSSION 

An interesting outcome of our analysis is the relative lack of sen- 
sitivity of the forecasted /nl constraints on the nuisance parame- 
ters: constraints on non-Gaussianity are tighter with larger statis- 
tics, rather than with a smaller samples with better controlled sys- 
tematics in cluster mass estimates. While this is the case for the 




0.795 



0.825 



Figure 10. Constraints at the 68 per cent confidence level on the /nl— 0"8 
plane for the bright subsamples of the three surveys, i.e. by including only 
clusters with fluxed 30 times larger than the detection flux. Dot-dashed 
cyan, dotted blue and solid red curves are for the Deep, Medium and Wide 
surveys, respectively, while the short-dashed green curve is for the combi- 
nation of the three surveys. We assume here strong prior on the nuisance 
parameters and we include FM from Planck experiment. 



purpose of constraining non-Gaussianity, it does not necessarily 
hold for constraints on the Dark Energy equation of state. In or- 
der to quantify the effect of an uncertain knowledge of the nui- 
sance parameters on Dark Energy constraints, we use as a refer- 
ence the redshift dependence of the equation of stat e, assumed e.g. 
in the Dark Energy Task Force (DETF) report (e.g. jAlbrecht et all 
2009), w(a) = wo + (1 - a)w a . Introducing the figure-of-merit 
FoMdetf = (det [Cov(wo, w„)])~ 1/2 , we find FoMdetf - 2066, from 
the combination of mass function and power spectrum, by adding 
the contribution of the three surveys and assuming strong prior on 
the values of the nuisance parameters. Such a high value drops to 
FoM D etf - 987 under the more realistic assumption of weak pri- 
ors on the nuisance parameters, while further reducing only by a 
small amount, FoMdetf - 972, if we assume no prior on such pa- 
rameters. Results on the constraining power of the high-sensitivity 
X-ray surveys reported here on different Dark Energy models will 
be presented in a future paper. 

A first word of caution in the interpretation of such high values 
of the DETF figure-of-merit lies in the assumption that scatter in the 
mass bias is always assumed to be log-normal distributed. As dis- 
cussed bv lShaw et alj (2010). deviations from this assumption are 
in general negligible as long as one relies on mass proxies which 
have a small value of the scatter, while the uncertain knowledge of 
the distribution of the scatter may significantly worsen constraints 
as the scatter increases. In our samples of clusters identified down 
to the detection limit, cluster selection is based on X-ray luminos- 
ity, whose correlation with cluster mass can easily have a scatter as 
large as 30 per cent. On the other hand, the effect should be much 
smaller for the bright samples, that only contains clusters for which 
measurements of robust low-scatter mass proxies are available. Re- 
stricting our analysis to the bright surveys, we find FoMdetf = 423 
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and 176 by assuming strong priors and no priors on the nuisance 
parameters, respectively. 

Another word of caution lies in the assumption that redshifts 
for all clusters in the surveys are perfectly known. An obvious ob- 
jection when discussing the cosmological exploitation of large sur- 
veys of clusters concerns the possibility of measuring redshifts for 
all of them. While discussing in detail the synergies between future 
X-ray and optical/near-IR surveys is beyond the aim of this paper, 
it is worth pointing out that future imaging and spectroscopic sur- 
veys, both ground based (e.g., BigBOSSjJ PANSTARlfl LSSlfl) 
and from space (e.g., EUCLIEQ, JDEIVo) will provide either pho- 
tometric or spectroscopic redshifts for all clusters identified within 
the WFXT surveys. Furthermore, all clusters in the bright sample 
will have a sufficiently large number of counts to guarantee a rather 
precise X-ray spectro scopic redshift meas urement from the posi- 
tions of metal lines dGiacconi et alj|2009l) . In general, uncertain- 
ties in redshif t measuremen ts could be taken into account in our 
analysis (e.g., ICunhdl2009T) . at least as long as one has a reliable 
estimation of redshift errors expected from different observational 
techniques. While redshifts from optical spectroscopy are gener- 
ally so accurate that any uncertainty can be neglected, the same is 
generally not true for photometric redshifts. We also note that no 
attempt has been pursued so far to assess in detail the reliability of 
measurements of clusters redshifts from X-ray spectroscopy, as a 
function of signal-to-noise, cluster temperature and redshift. 

In our analysis we did not include the effect of red- 
shift space distortions in the dis tribution of galaxy clusters in- 
duced by peculiar velocitie s (e.g.. iKaisejl 19871 ; IWhite et al]|2009l ; 
iDesiacques & Shetr]|2010l) . The study of this effect to available 
galaxy redshift surveys (e.g jGuzzo et alj2008l) has indeed demon- 
strated that it provides importa nt constraints on the growth rate of 
density perturbations (see also Linder 2008). Although we expect 
that its applications to cluster surveys may be limited by the sparser 
sampling offered by the cluster distribution, the effect of redshift- 
space distortions should be in principle included when forecasting 
the cosmological constraining power of future cluster surveys. We 
will present this analysis in a forthcoming paper. 

The analysis presented here demonstrates the potential that 
future high-sensitivity X-ray cluster surveys could have in con- 
straining possible deviations from Gaussianity. However, it may be 
worth asking what current data can tell us about such deviations. 
For instance, a positive skewness ha s the effec t of an ticipating the 
first collapse of massive DM halos. |jee et alj {2009) recently re- 
ported the discovery of an unexpectedly massive galaxy cluster at 
z =! 1.4, XMMU-J2235.3, identified as part of the initial 11 sq.deg . 
of the XMM Distant Cluster Project survey dMullis et all 12003) . 
having a flux limit o f 10 x 10~ 14 erg s^'cirT 2 . Base d on weak lens- 
ing dJee et al.ll2009l) and X-ray dRosati et al.ll2009l) analyses, a ro- 
bust lcr lower limit of 5 x 10 14 is obtained for the cluster virial 
mass. By assuming a WMAP-5 cosmology, with cr 8 - . 81 an d 
Q„, = 0.28, and using the mass function by I Jenkins etafl feOOlh . 
Jee et al. (2009) found that only ^ 5 x 10 3 of such massive clus- 
ters should be expected within the survey area. Thus, they con- 
cluded that XMMU-J2 235.3 is a rather unlikely event in a standard 
cosmological scenario. Ijimenez & Verde] d2009T) argued that, for a 



fixed value of erg (= 0.77 in their analysis) the expected number 
of such massive clusters can in fact be significantly enhanced in the 
case of a positively skewed non-Gaussian distribution of primordial 
perturbations. 

In Figure [JTJ we show the curves in the a 8-/nl plane corre- 
sponding to different numbers of clusters expected at z > 1 .4 within 
1 1 sq.deg. and having mass of at least 5 x 10 l4 M o . R esults are given 
for the reference non-Gaussian mass function from ILo Verde et all 
(2008) (dot-dashed cur ves), that we used for our forecasts, and 
for the mass function bv lMatarrese et al] d2000h (solid curves). For 
both mass function s, we applied the correction to A c suggested by 
iGrossi etal1d2009l) . As previously discussed, these two mass func- 
tions come from different approaches to approximate the exact re- 
sult for small values of /nl- As expected, the difference between 
the two mass functions becomes non negligible for / NL > 100 for 
the rare event of such a massive cluster at z — 1.4. For each model, 
the four curves, from right to left, are for 0.05, 0.02, 0.01 and 0.005 
such massive clusters found within the survey a rea, respectively . 



For homogeneity with the analysis carried o ut by I Jee et al. 
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we used here the Gaussian mass function by Ijenkins et al 
While / NL and <r 8 are left free to vary, all the other cosmological 
parameters are kept fixed at the fiducial values adopted in our ref- 
erence cosmological model (see Sect. [3)- The results shown in this 
plot confirm that a positive skewness helps increasing the expected 
number of clusters. The effect of non-Gaussianity is strongly de- 
generate with that of changing erg. For instance, increasing the ex- 
pected number of clusters by about a factor of ten for cr 8 = 0.8 re- 
quires / NL valu es in excess of the ran ge allowed already at present 
by CMB (e.g.. iKomatsu et ak 20 id and refer ences therein) and 
Large Scale Structure (LSS) " dSlosar et al.ll2008h data. On the other 
hand, the same boost in the cluster number can be achieved by re- 
quiring /nl - 100 and increasing erg to a 0.87, again in tension 
with current CMB and LSS constraints. The conclusion of this anal- 
ysis is that for XMMU-J2235.3 not to be a very unlikely event, a 
degree of non-Gaussianity in excess of the currently allowed CMB 
bounds is required, unless one wants to violate current constraints 
on o"8. Clearly, more than a single detection of such massive distant 
clusters are needed to draw firm conclusions. However, this exam- 
ple further confirms the strong constraining power of even few mas- 
sive clusters at z > 1. In addition, since galaxy clusters probe much 
smaller scales than the CMB, they offer a complementary approach 
to test a possible scale-dependence of non-Gaussianity. 



6 CONCLUSIONS 

We presented forecasts on the capability of future high-sensitivity 
X-ray surveys of galaxy clusters to provide constraints on de- 
viations from Gaussian primordial perturbations. Our analysis is 
based on computing the Fisher Matrix (FM) for the information 
given by the evolution of mass function and power spec trum of 
galaxy clusters. Following the approach by Tegmark (1997) t o com- 
pute the power-spectrum FM (see also iMaiumdar & Mohd 120041 : 
Rass at et al . 2008; IStril et ai]|201Ch . we include in the analysis the 
information related to the possible scale-dependence of the lin- 
ear bi as, which re presents a unique fingerprint of non-Gaussianity 
(e.g.,|Verde 2010, and references therein). According to the self- 
calibration approach, the model parameters entering in the FM esti- 
mate are nine cosmological parameters and 4 nuisance parameters, 
the latter defining the relation between cluster mass and observ- 
able upon which cluster selection is based. Our analysis is based on 
assuming an observational strategy designed for the Wide Field X- 
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Figure 11. The number of clusters with mass larger that 5 X 10 I4 M o , found 
in the redshift range 1.4 < z < 2 within the same surve y area of 11 sq.deg. 
where the XMMU-J2235.3 cluster has been detect ed free et alj2009h . Non- 
G aussian mass function b y iLo Verde et alj (2008) (dot-dashed curves) and 
by Matarrese et al. 12000) (solid curves) are shown. From right to left, ma- 
genta, blue, red and green curves show the models on the o"s-/nl plane 
predicting 0.05, 0.02, 0.01 and 0.005 clusters within the survey area, re- 
spectively. All other cosmological parameters have been kept fixed to our 
reference values (see text). 

ray Telescope (WFXT, e.g. lGiacconi et al] |2009), in which a Wide 
Survey covering most part of the extragalactic sky is complemented 
by a Medium and by a Deep Survey (see Table QJ. The latter pro- 
vides mass proxies down to the flux limit for cluster identification 
in the Wide Survey (see Table 1). We showed forecasts for the two 
parameters that, for a fixed expansion history, define the timing of 
cosmic structure formation, namely erg and /nl, while marginaliz- 
ing over all the remaining parameters. Informations on such con- 
straints are quantified by introducing the figure-of-merit for struc- 
ture formation timing of Eq.J24t. 

The main results obtained from our analysis can be summa- 
rized as follows. 

(a) Power spectrum and number counts of galaxy clusters are 
highly complementary in providing constraints: while the former 
is sensitive to deviations from Gaussianity, through the scale de- 
pendence of the bias, the latter is mostly sensitive to cr 8 . 

(b) Most of the constraining power for these two parameters lies 
in the Wide Survey, while the Medium and the Deep Surveys play an 
important role for the estimate of X-ray mass proxies for 2 x 10 4 
clusters out to z ~ 1.5. 

(c) Combining number counts and power spectrum information 
for the three surveys turns into A/nl - 10 for the lcr uncer- 
tainty with which a deviation from Gaussianity associated to a 
"local shape" model can be constrained. Correspondingly, we find 
FoMsft - 39 for the figure-of-merit of structure formation timing. 

(d) Quite interestingly, while the value of FoMsft significantly 
worsens when assuming more conservative priors on the nuisance 
parameters, the above constraint on / NL is weakly sensitive on such 
priors. 



(e) Th e presence of a cluster as massive as XMMU-J2235.3 at 
z - 1.4 dJee et alj2009h turns out to be a rather unlikely event, even 
allowing fo r an amount of non-G aussianity consistent with c urrent 
CMB (e.g jKomatsu et alfeOld) and LSS dSlosar et alj2008h con- 
straints. This further demonstrates the strong constraining power of 
detecting an even small number of massive high-z clusters. 

Our analysis lends support to the important role that fu- 
ture cluster surveys will play in constraining deviations from the 
Gaussian paradigm, with far reaching implications on the primor- 
dial mechanisms which seeded density inhomogeneities. The re- 
liability of our forecasts relies on the possibility of calibrating 
to high precision a universal expression for mass f unction and 
large-scale bias. A number of independent groups dGrossi et al.l 
20071: iDalal et al.ll2008l;lDesiacques et alJl2009l:lGrossi et alj|2009l : 
Giannantonio & Porcianil 120091 : IPillepich e t al.l 120101) carried out 
large N-body simulations with non-Gaussian initial conditions, 
finding in general a quite good agreement for the calibration of both 
mass function and bias. However, the precision required for the cal- 
ibration of such quantities, for them not to spoil the constraining 
power of large surveys, is probably higher than what reached at 
present. First assessments of the impact of uncertainties in the mass 
func tion calibration on DE constraints have been already presented 
(e.g. lWuetal1l201Ch and indicate that such uncertainties may not 
be negligible. There is no doubt that larger suites of non-Gaussian 
simulations are required to calibrate mass function and large-scale 
bias also for a range of models beyond the local non-Gaussian 
models that we considered in the analysis presente d here. 

A few days after our paper, ICunha et alJteOlOl) also subbitted 
a paper regarding the study of constraints on non-Gaussian param- 
eter f NL from cluster surveys. They use the Fisher matrix approach 
applied to the count-in-cell method to extract information on evolu- 
tion of cluster nu mber density and clustering. Differ ently from the 
previous paper bv lQguril ( l2009h . lCunhaetal.l bOlfj) also included 
the contribution from covariance between counts in different cells. 
This allowed them to sample the power spectrum ov er a large scal e 
range, thus obtaining tighter constraints on f NL than lOguril j2009h . 
By specialising their analysis for the mass selection and sky cover- 
age expected for the DESj optical survey, they forecast a precision 
of (Tf NL =: 1-5. A detailed comparison between our and their anal- 
ysis is not straightforward. Besides using different survey specifi- 
cations, our and their analyses also uses different prescriptions for 
the mass function and the bias. Just as an example of how sensitive 
the choice of the bias m odel is, we verified that excluding the bias 
correction suggested by iGrossi et al. | J2009l) (the factor q = 0.75 
that we introduce after Eqll U. the expected errors in f NL would de- 
crease by approximately a factor of 2. This further emphasizes the 
need for a precise calibration of the model mass function and bias 
through extended sets of non-Gaussian N-body simulations. 
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